Transcriptome profiles of metamorphosis in the ornamented pygmy frog Microhyla fissipes clarify the functions of thyroid hormone receptors in metamorphosis

Anuran metamorphosis is an excellent system in which to study postembryonic development. Studies on Xenopus (Mesobatrachia) show that thyroid hormone receptors (TRs) regulate metamorphosis in a ligand-dependent manner by coordinating the action of hundreds of genes. However, whether this mechanism is conserved among amphibians is still unknown. To understand the molecular mechanism of this universal phenomenon, we report the transcriptional profiles of the three key developmental stages in Microhyla fissipes (Neobatrachia): premetamorphosis (PM), metamorphic climax (MC) and completion of metamorphosis (CM). In total, 2,293 differentially expressed genes were identified from comparisons of transcriptomes, and these genes showed stage-specific expression patterns. Unexpectedly, we found that TRα was highly expressed in Xenopus laevis and Bufo gargarizans at premetamorphosis but showed low expression in M. fissipes. In contrast, TRβ was highly expressed during metamorphosis in M. fissipes and X. laevis. This result may imply that TRβ is essential for initiating metamorphosis, at least in M. fissipes. Thus, our work not only identifies genes that are likely to be involved in Neobatrachia metamorphosis but also clarifies the roles of unliganded TRα in regulating tadpole growth and timing of metamorphosis, which may be conserved in anurans, and the role of liganded TRβ in launching metamorphosis.